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Abstract 

We consider the linear response of a system modelled by continuous- 
time random walks (CTRW) to an external field pulse of rectangular 
shape. We calculate the corresponding response function explicitely 
and show that it exhibits aging, i.e. that it is not translationally in- 
variant in the time-domain. This result differs from that of systems 
which behave according to fractional Fokker-Planck equations. 

PACS numbers: 05.40.-a, 02.05.-r 

The response dynamics of several physical systems to external fields is 
well-described by the continuous-time random walk (CTRW) model intro- 
duced by MontroU and Weiss [Q. This model was extremely successful in 
the explanation of dispersive transport in amorphous semiconductors 0, 
see [0,^ for reviews. Recently it has been shown that in time- independent 
fields CTRW dynamics is sometimes rendered well through fractional calcu- 
lus relations such as fractional Fokker-Planck equations (FFFE) [Q. Now 
FFFEs have considerable mathematical advantages when compared to dis- 
crete stochastic schemes such as CTRW, since they allow to incorporate 
readily different initial and boundary conditions. In time-independent fields 
the FFFE expressions follow from those of CTRWs using the Kramers- Moyal 
expansion |^. Here we show that care is to be taken when treating time- 
dependent external fields. We demonstrate that for these the linear response 
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obtained in CTRW and in FFPE schemes may differ considerably; we display 
this here by taking exemplarily an external field which is switched on and off, 
i.e. which is of rectangular shape in time. In this case the CTRW-process 
displays aging, i.e. the response depends explicitly on the delay between the 
time of measurement and the time at which the system was prepared in a 
given state; hence the response reveals an essentially non-equilibrium pro- 
cess. On the other hand, FFPE follows the usual linear response behavior of 
systems close to equilibrium and is expected to be applicable in the cases 
of polymers [§] and rough interfaces . The differences between CTRW and 
FFPE turn out to be very pronounced in the time just after switching off the 
external field. 

As is usual in CTRW we envisage an ensemble of noninteracting particles 
which may be influenced by external fields (say, the particles are charged). 
The particles follow then CTRWs, i.e. sequences of jumps. The time intervals 
ti between the jumps are uncorrelated. Of interest are waiting times which 
follow power-law distributions, i.e. 

iPit) = 7/(1 + t/Toy+\ with < 7 < 1. (1) 

The physical motivation for such iplt) -ioims may be rationalized using ran- 
dom traps whose energy distribution is exponential |T^. In what follows we 
put To = 1 and work in dimensionless time units. 

A basic quantity in the CTRW formalism is Xn{t), the probability to make 
exactly n steps up to time t. In the standard, decoupled CTRW picture (in 
which the spatial transition probabilities between the lattice sites are inde- 
pendent of the waiting-times) the probability distribution P{r,t) of finding 
a particle at r at time t given that is started at at time obeys [ p!T[] 

00 

P(r,t)=^P„(r)x„(t), (2) 

where -Pn(r) is the probability to reach r from in n steps. We note that 
Eq. (^ is the starting point for the derivation of FFPE in Ref. . 

Using Eq.(§) several important relations follow: Thus the mean particle's 
displacement equals X(t) = J2'i^=o^nXn(t), where X„ is the mean displace- 
ment after n steps. In a weak homogeneous and time-independent external 
field E, X„ is proportional to E and to n. In particular, if only jumps be- 
tween nearest neighbor (nn) sites are allowed, the mean displacement per 
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step is 

^ = Er.^^ = AE, (3) 

i B 

where the sum runs over all nn vectors, is the Boltzmann constant, T is 
the temperature and fi is the mobility tensor. In highly symmetrical lattices 
X II E for any field direction and Eq.(^ reduces to a scalar relation, x = /iE. 
Furthermore, since X„ = nx =n/iE: 



X{t) = J2 ^nxn{t) = xN{t) = fiENit), (4) 

where N(t) = J2iZo f^Xnif) is the mean number of steps performed up to time 
t. 

In what follows we restrict ourselves to isotropic lattices and use a scalar 
notation. Typically, the evolution of N{t) depends now on the initial condi- 
tions 0. Thus, if we are interested in the number of steps performed from 
the moment t = in which the system was initially prepared (this is the 
usual situation in photoconductivity), all the intervals between the steps are 
governed by the same ip{t). Then A^(A), the Laplace-transform of N{t), reads 



A 



as is immediately evident from the expression for Xn(A) = ip{X)'^ [1 — V'(A)] / A, 
Ref. [|T^]. The Laplace transform of ip{t)^ Eq.(|I]), for small A is known to be 
il'{X) = 1 — A'^r(l — 7) ||12[. From Eq.(^ we now have 



m ^ Y^f'^^'' (6) 

The inverse Laplace-transform of this expression reads: 

Nit) = ^- -t^ = (7) 

^ ' r(l + 7)r(l-7) 7r7 ^ ^ 

Consider now the case that E is switched on at the time after the 
initial preparation of the system. In this case one has x = for t < t^, so 
that only the number of steps performed after tyj matters. This situation was 
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discussed in Ref. (see for a review). Now, the mean number of steps 
performed during the observation time r = t — t^,, i.e. between ty^ and t, is 

N{T,t^) = N{t)-N{t^) = '-^{t^-tl) = 

= "^[{r + t^y-tl]. (8) 

Eq.(^ fits very well the results of numerical simulations, see Figs.l and 2, 
where the simulated points are obtained from 10^ realizations of the random 
motion of walkers whose waiting-time distributions follow Eq. (|1|) . The results 
of the simulations are marked by symbols, and the results of Eq. (^ are shown 
as lines. The parameters of Fig.l are 7 = 0.5 and tw varies from 300 to 30000. 
In Fig. 2 we show the results for 7 = 0.3, 0.5 and 0.75 as functions of r for 
a fixed value of = 1000. Note that no fit parameters are used. 

Let us now return to the motion in the field E. The particles' displace- 
ments grows as fiENij) and the related current is 

J = f,E^N{T) = '^i,E{r + t^r~\ (9) 
dr IT 

Note that the overall response function a{T, t^) defined through j = cr(r, t,,u)E 
is given by cr(r, t.^,) = (vr"^ sin 717) fifl^^il + r/t^)'^''^, and shows simple scal- 
ing with respect to its two time variables: cr(r, t^) = t1^^F{t/t^). Thus 
in the limit of short observation times, t <^ t^,, the response function 
o'iTytyo) = (vr"^ sin 717) fit^^ is r-independent and hence describes an Ohmic 
transport. However, the value of the conductivity decays with the delay 
time. Systems in which the response to an external agent depends explicitly 
on the delay between preparation time and acting (measurement) time are 
referred to as aging systems. This kind of behavior was found to be very pro- 



nounced in CTRWs with < 7 < 1 ||14|-[T7[|. For t ^ tw, on the other hand 



j = (vr ^ sin 777) fiEr"' ^; the current is hence dispersive and independent of 

Let us now consider the response of the system to a field switched on at 
time tyj and switched off at time tz- Note that when the field is switched off, 
the directed component of motion ceases immediately: there is no afteraction. 
Hence the current is 

t <tw 

j(t) = { (tt-i sin7r7) fiEf^-^ t^<t<tz (10) 
t>t. 
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which is a causal response concentrated on the time interval in which the 
field acts. 

We now continue by discussing our findings in relation to a recently intro- 
duced approach to slow relaxation in time-independent external fields, based 
on fractional Fokker-Planck equations (FFPEs) P,|l9| or on fractional Mas- 
ter equations |2^. The FFPE describing subdiffusive behavior in an external 



field reads: 

^P(x,t)= oDl-''UpP{x,t), (11) 

where P{x,t) is the pdf to find a particle (walker) at point x at time t. 
In Eq. (pi]) oDI~^ is the fractional Riemann-Liouville operator (0 < 7 < 1) 
defined through [§] 

a-7...^_ 1 d f^^^, Z{t') 



* ^ ' r(7) dt Jo {t - t'y-^ ^ ' 

(in which t = can be associated with the time at which the system was 
prepared). Furthermore, in Eg . ([TT|) £i?p is the Fokker-Planck operator 

Cfp = - /i|-/(x) = - f,E^. (13) 

ax^ ox OX'' ox 

The r.h.s. of Eq. (|l^) holds, since in our case the acting force f{x) = —dU / dx 
is homogeneous, x-independent and of magnitude E. The fractional Fokker- 
Planck equation, Eq.(|ll]), has turned out to be useful in describing many 
phenomena connected with anomalous diffusion or relaxation patterns P]. 
We note that Eq.(|rT|) can be derived from the decoupled scheme, Eq.(^ 
using the Kramers-Moyal expansion, Ref. 0, provided the field is time- 
independent. The free relaxation properties of the CTRW system are re- 
produced by Eq.(|lT]) and can be expressed in terms of the Mittag-Leffler 

functions [|,|2D|. 

We can now make connection to X{t) of Eq.(|D by remarking that 



X{t) = / xP{x,t)dx. (14) 

Then multiplying Eq.([Tl|) by x and integrating it over the whole axis we get 
(parallel to a simple diffusion equation): 

d f°° i_ ( f°° d"^ f°° d ^ 

— / xP{x,t)dx = o-Dt ^IK x-—^P{x,t)dx — fiE{t) / x—P{x,t)dx 
ot J— 00 y J —00 ox J— 00 ox J 

(15) 
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The left hand side of Eq.(|r^) is nothing but ^X(t), whereas the right hand 
side can be simphfied by integration by parts. Using now that P{x, t) and 
its derivatives with respect to x vanish at infinity (as is reasonable for our 
models), the first integral vanishes and the second one is unity. Hence 

j^)= oDl-'fiEit). (16) 

Concentrating on the temporal dependence, we get in the case that the 
external field is switched on at t = and off at t = t^, i.e. E{t) = 
E9{t' - t^)0{t^ - t') that 

■^^ ' dt ^ ' r(7) dtJo (t-ry-^ ^ ' 

Evaluating the r.h.s. of Eq.(16) we get: 

t <tyj 

j{t) = \ [^l/T{^)]E{t-t^y''^ t^<t<t, (18) 



t > t,. 



which shows a functional dependence on the external field that strongly dif- 
fers from Eq. ([TO|) as soon as t > t^. 

Equations ( p^Of ) and (|1^) coincide only if one supposes /i/r(7) = jin^^ sin 777 
and takes = 0, — > oo\ this limit parallels the findings of Ref. 0. In 
general, however, the difference is large: Thus, Eq.(|l^) describes a response 
concentrated on the time- interval of the field- action, tw <t <tz'- no afterac- 
tion effects are seen. The current never changes sign and has finite jumps at 
t = and t = tz- On the other hand, Eq . (|18D shows considerable afterac- 
tion: the current does not vanish for t > t^. Moreover, the current diverges 
at t = ty^ and t = and changes its sign from positive to negative at t = t^. 
Moreover, the overall response described by Eq.([TE|) is invariant under time 
translation, i.e. depends only on the differences t — t^, and t — tz, which is 
not the case for CTRW, Eq.(|ig). 

The structure of the response, Eq.(0) derives from linear response close 
to equilibrium where in general the response I{t) and the field E{t) 

are related through the causal linear integral operator 



/(t) = f (j){t-f)E{t')dt' (19) 
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which does not account for aging effects. Memory effects in CTRW do not 
arise from a memory kernel, as in FFPE, but from temporal subordination 
(see Refs. |T0|, p!^] ) , which is an integral construct different from the convolu- 
tion of Eq.(^). Note that systems whose dynamics is governed by sub diffu- 
sive CTRW (i.e. for which 7 < 1) show aging; thus measuring the system's 
response to a pulsed field can be instrumental in determining whether the 
system obeys CTRW, Eq.(|ig) or FFPE, Eq.(|l8|). 

Let us summarize our findings. We have considered the linear response 
of systems governed by CTRW and by FFPE dynamics to an external field 
switched on at t = t^. We showed that for CTRW the response is initially 
ohmic and that it becomes dispersive after times comparable to the waiting 
time t^. We also studied the CTRW response to a rectangular field pulse 
of finite duration, switched on at t = ty, and off at t = t^. The form of 
this response (retardation effects after switching the field on, but absence of 
afteraction after switching it off) shows a behavior which differs considerably 
from that displayed by systems obeying FFPE. 
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Captions 



Fig.l. Mean number of steps N{t) performed in time t, displayed in 
double-logarithmic scales, for CTRW with a waiting-time distribution, Eq.(l) 
with 7 = 0.5. Here the triangles, circles, diamonds, squares and crosses 
indicate the t^-values 300, 1000, 3000, 10000 and 30000, respectively. The 
full hues reproduce Eq.(|). The dotted line has the slope 1, the dashed line 
has the slope 1/2. 

Fig. 2. Aging effects, to be seen from the dependence of on at 
t = 1000. The CTRWs waiting time distributions obey Eq.(l): we consider 
three different values of 7; namely 7 = 0.75 (circles), 7 = 0.5 (squares) and 
7 = 0.3 (triangles). Note the double-logarithmic scales. The full lines again 
reproduce Eq.(g). 
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